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Abstract 

A new method based on the Monte-Carlo calculation on the lattice is proposed 
to study the Casimir effect in the noncompact lattice QED. We have studied the 
standard Casimir problem with two parallel plane surfaces (mirrors) and oblique 
boundary conditions on those as a test of our method. Physically, this boundary 
conditions may appear in the problem of modelling of the thin material films interac- 
tion and are generated by additional Chern-Simons boundary term. This approach 



(N 

> 
O 



Oh! 



> 

in 



o 



for the boundary condition generation is very suitable for the lattice formulation of 
^ . the Casimir problem due to gauge invariance. 
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1 Introduction and general motivation 

During the last few years Casimir effect has attracted much attention due to the great exper- 
imental and theoretical progress in studying of this phenomenon. This macroscopic quantum 
effect plays crucial role in nanophysics, micro-mechanics, quantum optics, condensed matter 
physics, material science, and also it is very important for different models of boundary states 
in hadron physics, heavy-ion collisions, and cosmology. 

Nowadays, many theoretical methods for calculation of Casimir effect are proposed. Un- 
fortunately, most of these methods are based on the fixed boundary conditions or external po- 
tentials. Such a simplification can lead to problems in gauge invariance, renormalizability and 
locality. Moreover, there are some problems with different thermal corrections to Casimir force: 
different methods of calculations of the Casimir effect predict different corrections. Another 
part of the problem is calculation of Casimir force for complicated forms of boundary surfaces. 
Typically, various approximate methods (like the proximity force approximation method PQ, 
[2]) are used in the case of curved surfaces, but it is still unclear now whether those approx- 
imations are correct and for what tasks they can be applied. And finally, analytic methods 
for Casimir effect calculation are very complicated and strongly dependent on the shape of 
surfaces. Practically, Casimir effect was studied analytically only for cases of plane, spherical 
and cylindrical surface forms. More complicated tasks have not been studied very well by this 
moment but such cases typically appear in experiments. 

Based on arguments discussed above, it seems a very important task now to create a general 
method for calculation of Casimir effects which would work well for different shapes of boundary 
surfaces, at non-zero temperature and density and under other external factors. It means that 
such a method should be formulated very generally for working in different coupling regimes 
and different external conditions. And we think that direct lattice calculations in quantum field 
theory can meet all these requirements. 
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In our paper we consider the simplest Casimir problem with two parallel planes as a test of 
our method. The existence of the analytical answer for this problem is the additional motivation 
for such a choice of bound surfaces. This analytical answer assists us in fitting procedure for 
our numerical results. 

A crucial obstacle on the way of the realization of any Casimir problem on the lattice is the 
following. If the Casimir energy is a reaction of the vacuum on the presence of the boundary, 
what is "the boundary" in terms of lattice formalism? In other words, what is an observable 
quantity corresponding to such boundary? In fact, the answer is non-trivial and we devote the 
first part of the paper to this question. In the second part, we discuss the lattice algorithms 
and numerical results. And the last part is a conclusion. 

2 Chern-Simons boundary conditions and Casimir effect 

Casimir effect is a reaction of the vacuum on boundary condition. The spectrum of vacuum 
fluctuations depends on the boundary conditions. Changing of the boundary conditions leads 
to changing of the spectrum of vacuum fluctuations and so to generating of the corresponding 
Casimir force on the boundary. In the standard quantum field theory formalism, such changing 
of spectrum of vacuum fluctuations can be described, for example, by means of Green function 
method pQ. This approach is a very powerful tool for studying many essential Casimir tasks 
PQ. Unfortunately, the application of this analytical method to the case of more complicated 
shape of the boundary surfaces is not so easy due to calculation difficulties. Our aim is the 
creation of the numerical method for the Casimir effect calculation directly from the quantum 
field theory action. The lattice formalism looks very attractive for this role but manifestly we 
can not base in our approach on the separation of vacuum modes corresponding with boundary 
from the full spectrum of vacuum fluctuation. In lattice formalism we work in Euclidian space 
and deal with full spectrum of vacuum fluctuations and can not easily snatch out vacuum fluc- 
tuations corresponding to some boundary conditions. We need some very delicate approach for 
separation of vacuum fluctuation modes that preserve gauge invariance of our lattice formalism. 
Fortunately, such an approach to Casimir problem was proposed recently [HUE]- 

This approach is based on a very elegant idea coming from some unique properties of the 
Chern-Simons action in three dimensions [5j[6]. Let us consider electro-magnetic fields in 3+1 
dimensions with the Maxwell action and additional Chern-Simons action given on 3-dimensional 
integral on the boundary surface S: 

S=~J d A x F^ - ^d 3 se^n a A^x)F up (x), (1) 

where e a ^ v P is the Levi-Civita tensor and n a is the normal vector to the boundary surface S, A 
is a real parameter. 

Let us consider now the simplest form of boundary surface S, namely two parallel infi- 
nite planes placed at the distance R from each other. The Chern-Simons formulation of this 
canonical Casimir problem was studied analytically in series of works [TDl [EE] • We will use this 
analytical answer for the fitting of our numerical data. 

In the case of plane form of the boundary surface S the Chern-Simons action in (TjQ) has the 
following form: 

S cs = \j (6(x 3 ) - 5(x 3 - R))e 3 ^A^x)F up (x)d A x. (2) 

where, in our formulation of this Casimir problem, normal vectors to the planes are turned 
in the opposite directions. This choice of the normal vector orientation corresponds to our 
renormalization procedure based on the connection between open and closed Casimir problems. 
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Figure 1: /(A). 



If the parameter A is small, electro-magnetic fields obviously don't feel any boundary and 
are free. What happens if the parameter A becomes large and tends to infinity and fields 
dynamics on the boundary surface S is determined by Chern-Simons action? Let us consider 
the equation of motion obtained from the action (pQ): 

UA» + X(5(x 3 ) - 8(x 3 - R))e^ up A a d v A p = 0. (3) 

At A — > oo, it is easy to obtain from ([3]) a corresponding boundary conditions on the surface S: 

£|||s = 0, #„| 5 = 0, (4) 

where H n and En are normal and longitudinal components of magnetic and electric fields 
correspondingly. These conditions mean the nulling of the energy flux of the electromagnetic 
field through the surface. 

This phenomenon corresponds to the well-known property of the Chern-Simons theory in 3- 
dimension, namely to topological vortex (strings) generation [HE]. As it was shown by Witten 
[9], Chern-Simons theory, which is a quantum field theory in three dimensions, is exactly 
solvable by using nonperturbative methods and the topological vortices play a leading role in 
this process. It was shown that the partition function of this theory depends on the topology 
of the surface S and gauge group only, and tends to the zero if A tends to infinity. There are no 
propagation modes in this case and if the boundary surface S is closed, the gauge fields inside 
and outside this surface become separated from each other by the topological vortex fields on 
the surface. In the case of the finite value of A, we have a non-trivial interaction between 
topological vortices and electro-magnetic waves [8]. For our practical purposes it is enough to 
know that the parameter A plays the role of the effective regulator of the penetrability of our 
boundary surface for electro-magnetic wave. 

The analytical answer of Casimir energy per unit area for two planes is given by Green 
function method is the following [101 [11] : 



where function /(A) : liniA_»oo /(A) = 1 (Fig. 1) can be written as: 
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In our paper we will study the Casimir energy per unit area behavior at small A, which is 
the following: 

£ cas = -^ + 0(A 4 ) (6) 

3 Wilson "bag" and numerical lattice simulation of Maxwell- 
Chern-Simons theory 

3.1 General statements 

In our paper we use the four- dimensional hyper-cubical lattice in Euclidian space-time and the 
simplest form of the action for the noncompact lattice QED: 

s=§EEO*>. 

where the link and plaquette variables are defined as: 

Qi,n{%) =eaA li , 

A^9i iU (x) = Oi tU {x + ft) - 9i >v {x). 

Here a is a lattice step and the parameter (3 = 1/e 2 . Physical quantities are calculated in the 
lattice formalism by means of field configuration averaging, where the field configurations (the 
set of all link variables) are generated with the statistical weight e~ s . 

We have clarified in previous sections that additional Chern-Simons action describes Casimir 
effect. In order to find a lattice description of Casimir interaction between boundary surfaces 
let us consider Wilson loop, which describes the interaction of charged particles. Wilson loop 
can be written in QED as 

W c = e 9§ c A ^ = e^ J ^\ (7) 

The exponent in ([Zj) is the additional term to the action. This term describes the interaction 

of the field with the current J^(x) = g § S(x — £)<i£ M of charged particle. Configuration 

c 

averaging of Wilson loop (W(R,T)) (where R and T are dimensions of the loop) converges in 
Euclidian time in the limit T — ► 00 to: 

(W(R,T))^Ce- y ^ T , 

where V(R) is the energy of interaction between charged particles. The same method can be 
used for calculation of Casimir energy by means of Chern-Simons action. 



4 




Figure 2: Wilson bag for two plane surfaces. 



Analogously to the description of charged particles interaction by ID integral along Wilson 
loop, we will describe the Casimir interaction of surfaces by corresponding 3D integral. The first 
problem is that for stationary objects the action ([2]) is an integral from t = — oo to t — oo, so 
by analogy to Wilson loop, we should enclose the surface of the integration in t-direction. The 
integration surface for two planes is shown in Fig. 2. This closing procedure can be performed 
both for plane surfaces and for any curved surface in 3-dimensional space. As the result of this 
procedure, so-called Wilson bag [HI H] can be obtained. It can be written as 



jAisFpo-dSu 



where £ is closed 3-dimensional surface in 4-dimensional space-time. The final conclusion is 
that Wilson bag is (by analogy to Wilson loop) observable quantity which gives us Casimir 
energy of the objects, defined by the surface of integration. For two planes we will calculate 
the following object: 

W Bag (R,T) = e lXS ^ T \ (8) 



where 



S(R, T) = / dt 




o 

R 

+ J dz 

o 

And in the limit T — ► oo: 




dxdydz(5(z - R) - 5(z))e 3upa A l/ F p(7 + 



dxdydt(5(t — T) — 5{t))e iupa A v F pa . 



(W Bag (R,T))^Ce' E -^ T . 



(9) 



The second problem is to rewrite the Wilson bag in terms of lattice objects (links and 
plaquettes). The product 



A F 

J^u- 1 pa 



can be exactly constructed only in non-compact QED, because just in the non-compact QED 
the link variable 8^ u is a lattice analogue of A u and the plaquette variable Ptflv is a lattice 
analogue of F pv . 

There are two requirements for lattice representation of the Wilson bag: 
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Figure 3: "Corner": the basic element for the building of CS lattice action. 

1) The whole integral for Wilson bag should be the gauge invariant quantity. 

2) "Locality". It means that in lattice representation of the product A u F pa , A u and F pa 
should be given in the same point x. This requirement is non-trivial, because 9 PiP(T gives the 
value of F pa in the center of the plaquette, but 9 U gives the value of A v in the center of the 
link. These are different points. 

3.2 Gauge invariance 

A basic element for building of the CS lattice action will be a "corner" (Fig. 3). It is a product 
of the link variable and one of the neighboring plaquette variables. This product is a lattice 
analogue of the product A v F pa . 

First we will show how to construct the gauge invariant quantity from these basic elements. 
Let's consider one CS plane closed over the lattice period in all directions (V3). The "naive" 
discretization of this quantity in Euclidean space-time is follows: 

Scs = ^2 e ^pan p 9^ v (x)9^ pcj {x + dx pr7 ), (10) 
where the vector dx pa appoints to the plaquette that attached to the link 9i jV (x): 

dX pcr Cpfj^ dptrP ^ 'p(T@ ' } 
Qpai bpai Cpo 0, 1. 

There are 8 variants how to join a link with one of the neighboring plaquettes, two of them 
one can find in Fig. 4. 

Now let us consider gauge transformation of the "naive" lattice CS action (TlO|) .The gauge 
transformation of the link variable is follows: 

9i,u — > 0i,v + ®(x + u) - a(x), 

and the action ( flOl) transforms into the following expression due to the closing over the lattice 
period: 

Scs -> £,j,upan li 9i :V (x)9 PjPrT {x+dx pcT ) + y^ j e^ vpa n^a{x){9 PtPa {x+dx pa -p)-9 p ^{x+dx p(J )). 
XGV3 X6V3 

(11) 

The second term in (jTTi) appears due to the gauge transformation. To formulate the lattice 
CS action in gauge invariant form, this second term must be equal to zero. It is possible to 
show that this term is equal to zero if plaquettes in it are the sides of one 3d cube for all lattice 
sites x and normal vectors n p . We should choose a variant to join a link and a plaquette which 
has the property given above. One example of such a join is shown in Fig. 5. 
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Figure 4: Two types of "corners". 

This combination of "corners" is organized as follows: for every 3d cube the "corners" of 
different orientations that lie inside it (it means that "corners" are built from edges of this 
cube) have theirs link free ends being at one lattice site. This structure can be obtained for 
example by the following choice of the vector dx pa : 

dx pa — —p — a, 

This choice provides the nulling of the expression: 

>pa (x + dx pu — v) — 9 PtPa (x + dXpa-)) — 0, Vx, p. 

Thus the Chern-Simons action obtained from such a system of "corners" is the gauge in- 
variant. 

One can perform the similar consideration for any closed 3-dimensional surface. An example 
of such a closing for the "corner" with some orientation one can find in Fig. 6. Every "corner" 
should be a part of such a chain and every chain should be closed. 




Figure 5: Combination of "corners" with different orientations that leads to the gauge invariant 
lattice Chern-Simons action. 
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Figure 6: Closing for the "corner" of one orientation. 

3.3 Locality 

The next point of our consideration is the " locality" . As it is well known, the plaquette variable 
9 P ,pa gives us the value of F pa in the center point of the plaquette, while the link variable 
gives us the A p in the middle point of the link. On the other hand, the averages 

l/4(6i, u (x) + 6 hv (x + p) + 6 hv (x + a) + 9 hv (x + p + &)) 

and 

l/2{e p>pa (x) + 6 PtPa (x + f>)) 

give us the values of A v and F pa in the center of the 3-dimensional cube determined by the 3 
vectors p, v and a based on one site x. 

It's obvious that the product of these two averages is the same as the sum of all 8 variants 
of "corners", in other words we should use all possible kinds to join of a link and a plaquette 
inside one 3d cube. These 8 variants are shown on Fig. 7 for one possible orientation of a 
"corner". For another two possible orientations one has a similar set of variants. 

Let us consider now all 3 orientation of "corners". For each of them there is the set of 8 
variants of possible join a link with a plaquette. On other hand one can say that we have 8 
systems of "corners" (3 "corners" in one system). Each of these systems is similar to the system 
of "corners" shown in the Fig. 5 (links' free ends of all 3 "corners" are collected in one vertex of 
a cube) and forms the gauge invariant quantity after summation over the closed 3-dimensional 
surface. 

Finally, our lattice representation of CS action is the sum over all possible variants of the 
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Figure 7: The structure of integral for "Wilson bag". To obtain the gauge invariant quantity 
we sum inside its 3-dimensional volume over all the products of plaquette and neighboring link. 
All these variants for one possible orientation of link and plaquette inside one 3-dimensional 
cube are shown here. 

"corners" for three different orientations inside the close 3-dimensional surface: 

Scs = g/3 ^ e^p^AifcXfyi/fcH^ 
xev 3 

(12) 

The factor 1/8 in (11 21) comes from the averaging over the 8 variants of the link and plaquette 
joining. The expression 6i tV 9 PtPCT gives us a 3 e 2 A u F p(T in continuum limit and by means of the 
factor (3 = 1/e 2 one can eliminate e 2 from the action. 

3.4 Deformed lattice 

To improve the accuracy we have used the lattice deformed in one direction. This formalism 
is quite similar to one that was used in article [H] but instead of t-direction we deform one 
of spatial directions (for defmiteness z-direction) . Briefly this formalism can be described as 
follows: 

All link variables are connected with the field by the standard way: 

9i tlt = eaA^fi = 1,2,4 
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except links in z-direction: 

9 ij3 = aeaA 3 . 

Here a is the lattice step in all directions except z and aa is the lattice step in z-direction. The 
action of this theory can be written in the following way: 

x fj,<v;fj,,i>=£3 x v 

The additional Chern-Simons action is not changed on the a such deformed lattice. If n M (normal 
4d vector to the 3d surface of "Wilson bag") is oriented along the z-direction, deformed links 
don't participate in "corners" and nothing changes. If n M is oriented along one of non-deformed 
axes, one deformed link always participates in link or plaquette variable of every " corner" . For 
example: 

fy,30 P ,i2 = e 2 aa 3 A 3 F 12 , 
^A,i3 = e 2 aa 3 A 2 F 13 . 

But we don't need to divide by a, because this automatic presence of a in this expression simply 
corresponds to change of elementary lattice 3-dimensional volume in this part of "Wilson bag". 
Before deformation it was lattice 3d cube with volume a 3 . After deformation it's rectangular 
parallelepiped. Its edges have lengths a, a and aa, so its volume is aa 3 . 



4 Results of numerical calculations 

In our numerical simulations we calculate the Wilson bag quantity (E]) with the lattice repre- 
sentation of Scs given by the formula ( |T2l) . The Wilson bag surface is closed over the lattice 
period in x- and y-directions. E cas is extracted from the limit ([9]). 

First we study the dependence of E cas on R (distance between plates). In the Fig. 8 one 
can see E cas (R) dependence calculated on the lattice with dimensions 12 in x and y-directions, 
24 in z-direction (in order to decrease correlations over the lattice period in z-direction) and 
72 in t-direction (in order to eliminate temperature effects). We did't use deformation in this 
calculation. Parameter A is equal to 0.02 here. 

For this calculation we use (3 = 4, but indeed in our Casimir calculations nothing depends 
on f3. Casimir effect without radiational corrections doesn't depend on the electron charge in 
the continuous theory. On the lattice the independence appears due to the following reasons: 
1) there is no phase transition in the non-compact lattice QED and j3 here only plays the role 
of the scale parameter for a numerical value of link variables; 2) we eliminate this dependence 
from Scs due to the multiplication by (3 in the final expression (TT2]) . The absence of the 
dependence was also tested numerically. So (3 can be chosen in rather a large interval according 
to calculational convenience. 

It's well known in the continuous theory [T] that E cas should depend on R as C/R 3 . So 
we'll try to fit the point sequence obtained from the numerical simulation by the function 

E cas = Pi + ||. (13) 

Our aim is, of course, the coefficient P 2 . P\ is an analog of something like self energy of charged 
particles. Such self energy can also be obtained in Wilson loop calculation. The insertion in the 
Figure 8 presents the full dependence from R = 1 till R = 6. The main graphic is the increased 
segment from R = 2 till R = 6. The red line is the fitting that was done for points R = 2... 6 
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on the model function (JT3il . One can see that numerical points at lattice distances R = 2... 6 
lie on the fitting rather well. The point R = 1 differs from the normal C/R 3 dependence due 
to lattice artefacts at small distances. 

The values of coefficients Pi and P 2 are given in the Figure 8. Now let us consider the 
physical meaning of the coefficient Pi. A full Casimir energy of the interaction between two 
planes have been obtained in these calculations. This energy is expressed in units a -1 , where a 
is the lattice step. If we take into account that the area of the plane is equal to (aN) 2 , where 
N = 12 is the size of the lattice in x and y-directions, we can write the physical (dimensional) 
value of the Casimir energy density: 

1 P 2 1 P 2 N~ 2 

E cas . phys - ~^J^y ~ -Jfap' 

After comparison between this formula and ([5]) we can conclude that P 2 N~ 2 should be equal 
to coefficient at R~ 3 in (jSJ) or in (jSJ) (A is small in our calculations by reasons discussed later). 
And these values are rather close. P 2 N~ 2 = —4.44 x 10~ 6 ± 7 x 10 -8 and the theoretical value 
of this quantity for A = 0.02 is -5.066 x 10~ 6 . 

In order to improve the accuracy of calculation and obtain more points for the fitting 
procedure, we consider the lattice deformed in z-direction. As the planes are separated in z- 
direction, we'll be able to study E cas (R) dependence more carefully. In the Fig.9 one can see 
Eca S {R) obtained on the lattice with following parameters: 

• 12 normal steps (the length of each step is a) in x and y direction. 

• 72 normal steps in t-direction. 

• 120 small steps (the length of each one is aa, a = 0.2) in z-direction. 

So the whole volume and dimensions of this lattice are the same that the previous one. The 
only difference is smaller step in z-direction. 

The fitting procedure was realized again starting from distance R = 2a = 10 (era) (2 non- 
deformed steps is equal to 10 deformed steps). As in the previous case, numerical points are 
very close to the fitting function. There is one peculiarity which we want to mention. The 
energy is expressed in the units a -1 , because we didn't change the step in t-direction, but R 
is expressed now in the units aa. So if the E cas (R) dependence is really close to the C/R 3 
function, P 2 coefficient in fitting function (TT51) on the deformed lattice should be divided by 
a 3 to obtain the P 2 coefficient on the non-deformed lattice. And these numerical values of 
coefficients really satisfy this condition. 

The next step is the analysis of P 2 dependencies on the parameter A and on the size of 
the lattice. The first one is presented in the Fig. 10. In accordance with ([6]), the calculated 
coefficient at R~ 3 is proportional to A 2 for small A. Unfortunately, numerical errors increase 
rapidly for large A and in present calculations we are limited to small values of A. 

The dependence of P 2 on the lattice size is shown in Fig. 11. P 2 is proportional to N 2 , 
where N is the size of the lattice in x and y-direction. So the calculated Casimir energy of two 
planes is really proportional to theirs area. 

The last point is the consideration of the continuous limit of our calculations. As we 
mentioned above nothing depends on the parameter /3, due to Casimir energy independence of 
electron charge. So the only considerable stage in continuous limit is to go to the "large lattice" 
limit where size of the lattice N —>■ oo. There is one feature in this continuous limit procedure 
and we think that this moment should be emphasized. The lattice step a disappears from our 
considerations because the final result of our calculation is the dimensionless coefficient at R~ 3 
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in expression for Casimir energy per unit area. So we don't need to define the physical volume 
of the lattice and the value of a is not important for us. Because of arguments discussed above 
the continuous limit for P 2 coefficient is not sensitive to the limit a — > and we have only the 
limit iV — > oo. E cas (R) dependence in Fig. 8. can be an illustration of this continuous limit: 
the correct dependence E cas ~ R ^ reveals only at lattice distances R ^ 2. 

Conclusions 

In this paper we have proposed the numerical method for the Casimir energy calculation based 
on the lattice simulations of QED. We have combined three ideas: 

• The generation of the boundary conditions by means of the additional Chern-Simons 
boundary action and new lattice discretization of this additional action. 

• The lattice "Wilson bag" concept (the lattice presentation of the closed 3-dimensions 
surface in Euclidian 4-dimensions space). 

• The improvement of the numerical results by using deformed lattice in space direction. 

The combination of the first two ideas is in fact a lattice definition of the respective quantum 
observable for the Casimir energy and for the Casimir interaction between surfaces and last 
one is a very important numerical tool to archive the necessary accuracy. 

We have tested our method in the simplest case of the Casimir interaction between two 
plane surfaces and have achieved results close to the analytical ones for this problem. 

Acknowledgments 

We thank Prof. V.V. Nesterenko, Dr. I.G. Pirozhenko, Dr. M.A. Trusov and Dr. V.N. 
Marachevsky for interesting discussions about the Casimir problem and Andrew Zayakin for 
essential remarks about vortex dynamics of the Chern-Simons theory. The resources of MSU 
supercomputer center were used in our calculations. The work is partially supported by the 
Russian Federation President's Grant 195-2008-2. 

References 

[1] M. Bordag, U. Mohideen and V. M. Mostepanenko, Phys. Rept. 353 (2001) 1. 

[2] B.V. Deryagin, I.I. Abrikosova and E.M. Lifshitz, Quart. Rev. Chem. Soc. 10 (1958) 295; 
Uspekhi Fiz. Nauk 64 (1958) 493. 

[3] N. Seiberg, Phys. Lett. 148 (1984) 456. 

[4] M. Luscher, Phys. Lett. B 78 (1978) 465. 

[5] M. Bordag and D. V. Vassilevich, Phys. Lett. A 268 (2000) 75. 

[6] E. Elizalde and D. V. Vassilevich, Class. Quant. Grav. 16 (1999) 813. 

[7] M. Marino, Rev. Mod. Phys. 77 (2005) 675. 

[8] P. A. Horvathv. lirXiv: 0704 .32201 [hep-th]. 



12 



[9] E. Witten, Commun. Math. Phys. 121 (1989) 351. 
[10] V. N. Markov and Yu. M. Pis'mak, J. Phys. A 39 (2006) 6525. 



[11] V. N. Markov and Yu. M. Pis'mak, arXiv:hep-th/0505218 



[12] A. Hosaka and H. Toki, Phys. Rept. 277 (1996) 65. 

[13] I. Y. Malakhov and K. A. Sveshnikov, Theor. Math. Phys. 132 (2002) 1201. 
[14] G. Burgers and F. Karsch Nucl. Phys. B 304 (1988) 587. 



13 




Figure 8: E cas {R) at (3 = 4.0 and A = 0.02. 




Figure 9: E cas (R) dependence at (3 = 4.0 and A = 0.02 on deformed lattice. Step in z-direction 
is era, in other directions it is a, a = 0.2. 
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Figure 10: dependence. Lattice is the same that was used in Fig. 1. 
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Figure 11: P2{N) dependence. N is the size of the lattice in x and y-directions. Lattice without 
deformations, (3 = 4.0 and A = 0.02. 
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